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The method proposed by the present authors to deal analytically with the problem of Anderson 
localization via disorder [J.Phys.: Condens. Matter 14 (2002) 13777] is generalized for higher spa- 
tial dimensions D. In this way the generalized Lyapunov exponents for diagonal correlators of the 
wave function, (ipn,m), can be calculated analytically and exactly. This permits to determine the 
phase diagram of the system. For all dimensions D > 2 one finds intervals in the energy and the 
disorder where extended and localized states coexist: the metal-insulator transition should thus be 
interpreted as a first-order transition. The qualitative differences permit to group the systems into 
two classes: low-dimensional systems (2 < D < 3), where localized states are always exponentially 
localized and high-dimensional systems (D > D c — 4), where states with non-exponential local- 
ization are also formed. The value of the upper critical dimension is found to be Do — 6 for the 
Anderson localization problem; this value is also characteristic of a related problem - percolation. 

PACS numbers: 72.15.Rn, 71.30.+h 



I. INTRODUCTION 
A. Experiment and theory 

Disorder leads to important physical effects which are of quantum mechanical origin. This has been revealed by 
Anderson pj in the study of a disordered tight-binding model. This problem has attracted great attention over many 
decades. A breakthrough came with the scaling theory of localization |2(. All states in a one-dimensional system 
(1-D) are localized, whereas in 3-D a metal-insulator transition occurs. An analytic solution is only known for the 1-D 
problem 0. Although there was no general analytical solution available, there was consensus that in 2-D all states 
are localized. 

For quite some time after the advent of the scaling theory, many believed it to be essentially under control. This view 
is less secure nowadays, in part because recent experiments have challenged conventional wisdom about disordered 
2-D systems. The 2-D case still presents a problem which has become apparent by experiments 0,0,E10- These 
experiments are still being discussed controversially. Experimental reality is certainly more complex than a simple 
tight-binding model, but these results provide a good reason for revisiting this classic theoretical problem. 

Recently we have been able to solve the 2-D case analytically 0- We have shown that in principle there is the 
possibility that the phase of delocalized states exists for a non-interacting electron system. For energies and disorder, 
where extended states may exist we find a coexistence of these localized and extended states. Thus the Anderson 
metal-insulator transition exists and should be regarded as a first order phase transition. Consequently we have 
returned to the old idea of Mott 0, ^| that the metal-insulator transition is discontinuous. This alternative idea 
was in its history completely abandoned with the advent of the scaling theory of localization. Thus we expect a 
renaissance of it. 

This result implies the failure of the scaling theory of Anderson localization. Although this paper is published @ 
and constitutes the basis for a new analytical investigation of the Anderson problem in the present paper for the 
case of higher dimensions (N-D problem) one has to accept the following: (i) the paper Q requires an independent 
confirmation, which requires a certain time; (ii) in the history of the problem one has developed quite a few conceptions 
and this leads to a critical attitude towards the new theory - i.e. there is a resistance - as it departs from conventional 
wisdom. This asks for a critical evaluation of the new theory and its results 0, 0] . One should acknowledge that the 
problem is many-sided and quite complex. Thus logical errors are possible which do not lie at the surface. Perhaps 
a reference to B. Pascal is appropriate here: A truth is so delicate that any small deviation from it leads you to a 
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mistake, but this mistake is also so delicate that after a small retreat you find yourselves in a truth again. Without 
an exact analytic solution a discussion cannot lead to firm results. 

Let us add at this point some reflections of a historical nature. It is well known that exact results in the field of 
phase transitions are exceedingly rare (so to say they are exceptions) and they exist only under certain conditions and 
limitations. The situation in the field of Anderson localization shows in this respect a certain similarity to the history 
of the theory of the second order phase transitions and of the Ising model. The phcnomenological theoretical work by 
Landau |13| and the numerical work by Bethe |14| were prior to Onsager's analytical theory |15| and they were well 
accepted. But the results of the analytical theory departed strongly from the phenomenological and numerical results 
in that it gave different critical exponents. But there are limitations to the analytical theory. The exact solution of 
Onsager for the Ising model is only possible for 2-D, without magnetic field and only for the case of an interaction 
between nearest neighbours. This is also the case in our present investigation. An exact solution is only possible (a) 
for the tight-binding approximation with (b) diagonal disorder, where (c) on-site potentials arc independently and 
identically distributed, i.e. for the conventional Anderson model. 

A contradiction between results from phenomenological theories (Landau, scaling theory) and analytical results 
(Onsager, the present work) have (Onsager) or might (the present case) lead to scientific advances. The metal- 
insulator transition is a phase transition. In the case of such phase-transitions the results tend to be unstable and 
depend on the approximations. An exact solution thus rarely shows an agreement with results from a phenomenological 
theory and numerical work, because the two latter methods make use of additional approximations and assumptions, 
whose consequences are hardly controllable. To return to the advances made possible by Onsager's analytical theory: 
V. Ginzburg got in the year 2003 the nobel prize for the work (together with L. Landau [16j)j in which he used 
the phenomenological theory of second order phase transitions. This all means simply that phenomenological and 
numerical results frequently require certain corrections, which is only then possible, when an exact solution is known. 

The main aim of the present article is thus the generalization of the mathematical tools of the previous article Q 
to the case of higher dimensional spaces and a physical interpretation of the new results. 

B. Structure of the present article 

The outline of the present article is a follows. The article opens the possibility to determine the phase diagram 
of the multi-dimensional Anderson localization analytically. This is based on the analytic calculation of Lyapunov 
exponents. The argumentation is given at the beginning of chapter 2 in detail. There we also give a short derivation 
of the equations for the N-D Anderson localization problem, and a necessary summary of the results of the article ||. 
The connection between the Anderson localization problem and signal theory is discussed and the most important 
concept of the proposed method - the filter H(z) - is defined. To understand fully these aspects of the theory a 
knowledge of the first paper § is recommended. The filter H(z) is generalized for higher dimensions D of the space. 
The investigation of its properties and the corresponding physical interpretation constitute the content of several 
chapters. The theory for high-dimensional systems (D > 4) is presented in chapter 3, whereas the theory for low- 
dimensional systems is given in chapter 4. An appendix deals with the mathematical conditions for the physical 
interpretation of the filter H(z) and a more detailed discussion of related aspects of the problem. Because of the 
length of the present article we omit both the mathematical details of the theory as well as an extended discussion of 
the implications of the results. This will be given in a future publication, where also a different procedure for part of 
the derivation is used instead of an average over initial conditions as done here. 

II. NEW METHODS 

A. Anderson localization and generalized diffusion 

Recently we have been able to solve the 2-D case analytically ||. The tight-binding equation in 2-D is solved for 
the wave function ip ntm and the second moments (diagonal correlators ("0„ jm )): 

Ipn+l.m (-^ £n.rn)^Pn,m — l,m ^ ^Pn,m+m' ; (1) 

m' 

where the summation over ml runs over the nearest neighbours of site (n, m) in layer n that is in a space of dimension 
p = D — 1 = 1. We assume taking the limit to infinite size L — > oo in p-space. The equation is solved with an initial 
condition 



1p0,m = 0, V>l,m = a rn- 



(2) 
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The on-site potentials £ n ,m are independently and identically distributed with existing first two moments, (£ n , m ) = 
and (e 2 n m ) =a 2 . 

The mathematical idea shortly discussed above is based on a simple physical idea. The aim of the article [2j 
consisted in determining a phase diagram for the Anderson localization problem. To this purpose the concept of the 
Lyapunov exponent 7 |8( is extremely well suited. The Lyapunov exponent plays the role of an order parameter: 
7 = for the metallic phase and 7^0 for the insulating phase. The value 7^0 for the insulating phase implies a 
divergence of certain mean values as a function of the index n. If we understand what this divergence means it will 
become clear, which quantities have to be calculated. 

This divergence is a general property of a number of different stochastic systems. To give a clear picture we consider 
in the following the equations for 1-D systems and in particular for diffusion (random walks) 

ip n+ i =ip n + £n, (3) 

and Anderson localization 



1pn+l = (E ~ E„)lp n ~ 1pn-l- (4) 

For symmetric diffusion, (e n ) = 0, (e^) = c 2 , the simplest and most important characteristic is the second moment 
(4>n) = V'o + & 2 n. A divergence of this variable for n — » oo is an unambiguous proof for the existence of diffusion. It is 
also analytically possible to calculate other even moments, e.g. (tpn)- These, however, do not supply new information. 
All even moments diverge simultaneously, and they do this as a power of n. A similar divergence also occurs for other 
mean values of the functions of \ip n \, e.g. (IV'nl)j however, these can only be calculated numerically. In other words, for 
the determination of the existence of diffusion it is fully sufficient to calculate only one moment, the second moment. 

The equation (J1J is nothing else but a generalization of equation J2J- The mean value (ipn) = fin) describes in 
this case generalized diffusion. If f(n) is bounded, f(n) < oo, then the proper dynamics in equation Q is stable 
and we have no diffusion (Lyapunov exponent 7 = 0). This corresponds in a physical interpretation to the existence 
of delocalized states. A divergence of the function f(n) for n — * 00 corresponds to generalized diffusion (localized 
states); in this case one could distinguish in addition between non-exponential localization (f(n) is a non-exponential 
function) and exponential localization (/(n) oc exp(27n) with 7 7^ 0). 

It is quite important to stress here that the formation of diffusion in a dynamic system represents always a type of 
phase transition. Therefore one can easily find the corresponding critical values and determine the phase diagram of 
the system. As for proper diffusion, an analytic calculation for eq. J3J is only possible for even moments [3, 8]. 

For these reasons it is fully sufficient to consider only the equations for the second moments (correlators) to determine 
the phase diagram of the system @ . This idea has a certain similarity to the method of Pendry 01 • Pendry has 
introduced the density matrix to the problem of waves in disordered systems. However, he was of the opinion (without 
giving a proof), that for Anderson localization all higher-order density matrices and hence all higher moments need 
to be considered. This program is practically unfeasible, so Pendry limited himself in to the second moment. 
Instead of solving the equations for the second moments exactly Pendry has mainly studied the properties of the 
transfer matrix and the process of its diagonalization. For the transfer matrix, however, the investigation of the 
thermodynamic limit is impossible (see discussion in 12]), and consequently in this way the phase diagram of the 
system cannot be obtained. 



B. Equations for correlators 



In the present paper we present the analytic solution for the general case D > 2. The knowledge of paper [8| is 
prerequisite for understanding the present one as it contains the full formalism. The generalization of this formalism 
to the N-D case presents no problem (see below). The main eqs. remain valid, only scalar quantities become 

vector quantities. 

The tight-binding equation in 1 + p dimension is (primed indices are summed) 

1pn-\-l,m ^n.mV^m ^Pn— l,m ~t~ ^~-m,m'^Pn,m': (5) 

ES m,m' / ^m+m" ,m' 7 (6) 

(summation over m" runs over the nearest neighbours) with initial condition *0o,m = and ^1 m — ct m . 
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One introduces the correlators (the averages are taken over disorder) 

U{n) m ,\ = (lpn, m 1pn-l,\) ■ (8) 

From the eq.® one easily obtains the system of equations: 

x(n + l) m .i = £ m ,icr 2 x(n) m! i + x(n - l) m j + (9) 
£ m , m 'a;(n) m ',i'£i',i _ An,m'l/(n)m',i - Ci,vy(n)v lJa , 

y(n + l) m .l = -J/(n)l, m + £m,m'ic(«)m',l- (10) 

They can be solved explicitly by introducing of the Z-transform 18J 

n*w = E^. (12) 

n=0 

which turn the equations into 

(^ - - V 2 8 m ,\)X(z) mil - l(l)m,l = £m,m'^^)m',l'A',l (13) 

->C m , m 'F(z) m ' ! i - C\,\>Y(z) v . m , 
zY(z) mA = -Y(z) hm + £ m , m >X(z) m >, h (14) 

Iteration of the second equation yields 

z 1 

Y(z) m ,\ = -n r£ m , m /.X"(,z) m ' 1 5 r-X"(z) m ,l'>Ci',l- (15) 

z A — 1 z z — 1 

We then obtain an equation for X(z) m j only, which can be solved by double Fourier expansion. A complete presen- 
tation of the method is rather lengthy and will thus be the subject of a following paper. An shortcut which quickly 
leads to the correct equations is suggested in section 3.2 of @): averaging over translations in p-space of boundary 
conditions, a m . 

Upon averaging over an ensemble of initial conditions (as described in |8|]) such that a m a m > = r m _ m / , the problem 
is translation-invariant in transverse directions. We then put: 

*(*)m,i = J -f±-X(z,k)e^™-V. (16) 

After averaging over initial conditions the diagonal correlator becomes independent of m: 

*^(^)m,m = ^ni (1*0 

f d p k 

X(z) m , m = X(z) = J ——X(z,k). 



(18) 



We obtain the final equations 



(Z ^'w 2 -£ 2 (k)}X(z,k) = T(k)+a 2 X(z), (19) 

p 

S(k)=E-2^2cos(kj), (20) 

- 2 = (21) 



(z + l) 



or 



X (o) (z)= (i±l) 



X(z) = H(z)X(°\z), (22) 



dPk r(k) 



(z-l)J (2tt)p [w 2 -E 2 {k)\ 
2 0+l) f d p k 1 



(23) 



H{z) (z-l) J (2tt)p [w 2 -£ 2 {k)] 



1 - a 7Z TT ,o.w-.2 c2.,.m . ( 24 ) 
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where the X^\z) (or xffl) refer to the ordered system (er = 0) and the X(z) (or x n ) to the disordered one (er =/= 0). 
For eq. ll'L'li the inverse Z-transform gives the convolution property 01 : 

n 

x n = ^atj ftn-j, (25) 

with 

oo , 
n=0 



C. Anderson localization and signal theory 



The essential poi nt in the analysis with respect to the localized or extended character of the states is to make use 
of signal theory |l8l | from electrical engineering and switch from an investigation of the moments x n to an analysis of 
the filter functions h n . 

In the theory of signals ^8|, x„ (or X^°\z)) is the input signal, and x n (or X{z)) is the output. Asymptotic 
behaviour of the solution is completely determined by the filter h n (or H(z)). The concept of the system function 
is a general and abstract description of the problem of localization. Thus the filter function has to be analysed to 
obtain general results and not the multitude of signals. This has been done in [8| for the 2-D case. This procedure 
has a certain similarity with the transition to an operator formalism in quantum mechanics. Particular signals x n 
depend on the initial conditions used and do not carry much physical information because of the unconventional 
normalization yj. For the localization problem the only property that matters is whether a signal belongs to the 
bounded or unbounded class and this can be derived from the filter. 

The essence of localization is contained in the filter H(z). We study the filter H(z) with properties described by 
generalized Lyapunov exponents 0, Il9| . The filter is a fundamental function of the disorder a only. 

A filter h n is uniquely characterized by a pole-zero diagram of its image H(z) which is a plot of the locations of the 
poles Xi and zeros of H (z) in the complex- z plane. We provide just a brief summary here, for more details consult 
HG3. The signals x { n ] and real, therefore H(z) will have poles and zeros that are either on the real axis, 

or come in conjugate pairs. For the inverse Z-transform H(z) h n one needs to know the region of convergence 
(ROC). Physical considerations dictate that only causal filters (h n = for n < 0) should be considered. They have 
ROCs outside a circle that intersects the pole with max |Aj|. A causal filter is stable (bounded input yields a bounded 
output) if the unit circle \z\ = 1 is in the ROC. Note that the explicit calculation of h n by the inverse Z-transform 
is not necessary, and it is also not feasible analytically due to the complexity of the function H(z). Only the type 
of the filter - stable or unstable - needs to be determined. The delocalized states (bounded output) are obtained by 
transforming the physical solutions inside the band (bounded input) provided that the filter H(z) is stable. Seeking 
for poles is quite a simple analytical task which gives rather general results by elementary methods. 

As an example for the general and abstract description of the problem of localization as stated above let us consider 
the following problem. It is well-known that disorders broadens the band; new states outside the old band arise for 
\E\ > Ef, — 2D. Are among these new states also extended states? Numerical work for D = 3 |20t l2l| ascertains this. 
This is the socalled reentrant behaviour of the mobility edge: the change from localized to extended states and back 
to localized ones upon increasing the disorder occurs for certain fixed energies. Because in the literature there is no 
physical explanation for this phenomenon one has simply accepted these results without critically examining them. 

With the help of signal theory we have found a particular transformation, which gives a connection between the 
states in an ideal system (zero disorder) and states in a system with disorder. Extended states (zero disorder) 
as input (bounded) signal transform into localized states (nonzero disorder) as output (unbounded) signal if the 
filter which is responsible for this transformation is unstable. If the filter is stable, then extended states (bounded 
input signal) transform into extended states (bounded output signal). It is known that for zero disorder outside the 
band, \E\ > Eb, there do exist only mathematical solutions which cannot be normalized. These correspond to an 
unbounded input signal. It is impossible to find a filter which permits a transformation of the unbounded input signal 
(mathematical solution) into a bounded output signal (extended states in \E\ > Eb). The reverse - the transformation 
of an unbounded input signal (mathematical solution) into a unbounded output signal (localized states in \E\ > Eb) 
- is on the other hand possible. I.e. the mathematical procedure developed by us generates in this case a negative 
answer to the posed question. 

We have shown in |8| that the filter H(z) is a non-analytic function of the complex variable z; this result remains 
valid also in the multi-dimensional case. The unit circle \z\ = 1 divides the complex plane into two analytic domains: 
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the interior and exterior of the unit circle. The inverse Z-transform is quite generally defined via contour integrals in 
the complex plane 

K = ±-ln{z)z n ^. (27) 

and this definition is only possible in an analytic domain. In this way in the formal analysis of the problem multiple 
solutions result. The first solution H + (z) is defined outside the unit circle and always exists. The filter H + (z) 
describes localized states and it is possible to connect its properties with the notion of the localization length 
The second solution H_(z) is defined inside the unit circle and does not always represent a solution which can be 
physically interpreted (this is the mathematical consequence that the filter be causal). The filter H_(z) describes 
delocalized states. The coexistence of the two solutions was physically interpreted in 0,0 as the coexistence of two 
phases - an insulating and a metallic one. Then the metal-insulator transition should be looked at from the basis of 
first-order phase transition theory. 



D. Conformal mapping 

The p-dimensional integral on the r.h.s. of eq. H24|l can be reduced to a one-dimensional integral. Consider the 
identity 



f 



£ 2 (k) " ;_«, w 2-(E + y y ■ (28) 

The integral representation of the Dirac ^-function and the Bessel function 

J W = [ e lxcos{k) dk. (29) 

will be used. Following let us define a complex parameter w — u + iv from eq. I|21|) in the upper half-plane, v > 0. 
Using the methods of complex variable theory we get 

' ' °' k 1 Y D (w,E), (30) 



(2tt) p J w 2 -£ 2 (k) iw 
where 

/•OO 

Y D (w,E) = [Jo^t)} - 1 cos(Et) exp(iwt) dt. (31) 
Jo 

Changing the complex variable z to the parameter w corresponds to the conformal mapping of the inner part [\z\ < 1, 
w = — [z 1 / 2 + z~ x / 2 ))] or the outer part [|z| > 1, w = (z 1 / 2 + z -1 / 2 )] of the circle onto the upper half-plane. The 
circle itself maps onto the interval [—2,2]. Note also that if H (z) = has complex conjugate poles, then on the 
upper w half-plane they differ only by the sign of u. To avoid complicated notations, we seek for poles in the sector 
u > 0, v > and double their number if we find any. 
The inverse function Q 

w 2 w i — 

z = -l + — ±-y/w 2 -4 (32) 

is double- valued. Its branch with the minus sign maps the w sector onto the inner part of the half-circle (|z| < 1). 
The second branch with the plus sign gives a mapping onto the half-plane with the half-circle excluded (\z\ > 1). 
Therefore in the parametric w-representation 

1 1 — A 



and 



(34) 
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III. HIGH-DIMENSIONAL SYSTEMS 



A. Upper critical dimension 

It is generally assumed that 2-D systems mark the borderline between high and low dimension ■ The existence 
of a transition in 3-D is not questioned (high-dimensional systems). These assumptions, however, originate from the 
scaling theory of localization. Here the marginal dimension is Dm — 2, and a phase transition exists only for D > Dm- 
Thus perturbation theory H for D = 2 + e (e <C 1) is possible. The effect of statistical fluctuations cause a change 
of regime at D c = 4 [1(1 . l23j |: in this way the upper critical dimension for localization D c arises. For D > D c there 
should not exist a phase transition. 

On the one hand these statements referring to higher dimensions are numerically nearly impossible to ascertain 
[24| . Statistics is bad and the length of the system L very small. Even for D = 3 progress towards extracting reliable 
numerical estimates of critical quantities has been remarkably difficult [2q|. On the other hand, if the results from 
the scaling theory of localization for 2-D systems are faulty (this is what we claim) , then the corresponding division 
of the systems into low and high dimensional ones is also wrong. Here one must develop an alternative picture. 

In the theory of critical phenomena |26ll27| many systems belong to a class, where an upper critical dimension D 
has a totally different physical meaning. It denotes the dimension, from where on the mean field approximation is 
exact, or where with other words all critical exponents reach stationary values. I.e. in this case the phase transition 
does exist also for D > D ; only in the limit D — > 00 the transition disappears, simply because the corresponding 
critical values have gone to infinity. Systems of quite different physical nature may have the same value of the upper 
critical dimension Dq. E.g. one finds Dq = 4 not only for the theory of magnetism |26J,|27j (in this case there exists 
also the marginal dimension, below which no phase transition is possible), but also in kinetics (cooperative phenomena 
in bimolccular processes by diffusion-controlled reactions) (2^, H^; in the latter case, however, there is no marginal 
dimension. Another example is percolation, where the upper critical dimension is Dq — 6 [30| . 



B. Stability and poles: solution H+(z) 

Let us start first from a purely mathematical comment: the integral Yd{w,E) is always finite for all w = u + iv 
in the sector u > 0, v > for high-dimensional systems with D > 4. This fact clearly follows from the asymptotic 
behaviour of the Bessel function for large values of its argument, Jo(2i) « — i= cos(2i — 7r/4). We shall further on see 
that the dimensionality D = D c = 4 is critical for localization, although this is no proper upper critical dimension 
D = D , which we shall determine below. Let us consider therefore the properties of the filter- functions in this region 
of the value of the dimension of space. 

Let us consider first the solution H + (z). According to the ROC of the causal filter is defined by the inequality 
\z\ > max I Ai|, where A^ are the poles. The case when the system function has a pole at z = A > 1 is the simplest one 
for an interpretation. In terms of signal theory [li§| the filter H + (z) is unstable since the pole lies outside the unit 
circle \z\ = 1 in the complex z-plane. As shown in [8j, such a filter describes exponentially localized states (insulating 
phase). In order to see this one can exploit a basic inverse Z-transform [is| : 

H(z) = z/(z - A) h n = A". (35) 

Therefore the pole of H + (z) at z = A = cxp(27) leads to exponential growth of the system function h n which in turn 
implies [eq. Q25JI ] an exponentially increasing mean squared amplitude x n . The growth exponent 7 is the so-called 
generalized Lyapunov exponent |3J related to the localization length by £ = 7" 1 . 

The value of the Lyapunov exponent 7 defines the phase. We start from the mathematical definition that all 
states with 7 7^ belong to the insulating phase. The states with 7 = on the other hand correspond to a metal. 
According to this definition the states with non-exponential localization also belong to the metallic phase, because 
they correspond to the value 7 = 0. We can consider these states as a bad metal, in contrast to a good metal, where 
one has truly delocalized states. 

We would like to give here a summary of the results which emerge from an analysis of the pole diagram (for details 
see the Appendix). The function Yd(w, E) defined by eq. (|31|l is purely imaginary for v — 0, u > uq, 

u =2p+\E\, (36) 

and the system function H + (z) itself is real. The pole must be located (if present at all) exactly in this region of the 
parameter w. It can be found as a solution of an in general transcendental equation 



a 2 n+{u) = 1, 



(37) 
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where 

O , ( u ) = — / [J (2t)] D ^ 1 cos(Et)sm(ut)dt. (38) 

Va 2 - 4 Jo 

For high-dimensional systems with D > 4 the function f2+ (it) has a maximum at u = uq and decreases monotonically 
for u > uq. The system function has a pole if the disorder exceeds a critical value, a > <jq(E), where 

a (E) =n+{u y 1/2 . (39) 

Therefore, in high-dimensional systems exponential localization takes place only if the disorder is strong enough. 

For a < <tq(E) the function H + (z) has no poles, its ROC its \z\ > 1. This means that the unit circle \z\ = 1 belongs 
to the ROC, and the filter is stable. 

We interpret solutions for this range of the disorder values, a, also as localized, however with non-exponential 
localization. Non-exponential localization corresponds to a Lyapunov exponent 7 = 0. This simply means that the 
corresponding function grows slower than the exponent exp(27n). To this class of functions do not only belong 
power-law functions, but also exponents with a different argument, e.g. exp(c • n s ) with 8 < 1. Here we encounter 
limits of applicability of the method, which come into play, however, only in physically inaccessible systems of high 
dimensionality. The elementary pole search can be applied only for exponentially localized states, the general case 
requires a detailed investigation of the filter. In its present form our theory is not capable to determine the type of 
non-exponential localization. 

The curve for <tq(E) (Fig. la) is the well-known mobility edge. For a > ctq{E) there exists an insulating phase, 
whereas for a < uq(E) a metallic phase is found (bad metal). A coexistence of phases is here not possible. Thus 
the phase transition here has an appearance as if it were a transition of second order. This simple idea, however, is 
contradicted by the behaviour of the Lyapunov exponent: the transition from 7 = to 7 7^ is not continuous. One 
can clearly see that all these mobility edges for high-dimensional systems show the same qualitative behaviour. 



C. Stability and poles: solution H-(z) 

Now let us turn to the second branch of the solution, H-(z). In this case existence of the poles leads to principally 
different consequences. Let us assume that the corresponding value of the parameter w is found and the pole z = Ai 
is located inside the unit circle, |Ai| = 1/A with A > 1. Formally, however, from the definition w — — (z 1 / 2 + z -1 / 2 ) 
the same value of w can be obtained for z = X2 = 1/Ai = A. The complex number A2 lies outside the region of 
definition of the solution, \z\ < 1. In this sense the pole at z = A2 is virtual. For the inverse Z-transform this fact is, 
however, irrelevant. The ROC for a causal filter is defined by the inequality \z\ > max|Ai| or \z\ > A > 1. Since the 
ROC and the region of definition of the solution \z\ < 1 do not intersect, a physical solution is absent. Therefore, the 
filter -ff-(z) as a physical solution is acceptable only if either there are no poles or they lie on the unit circle. The 
latter case is realized for D = 2 || and corresponds to so called marginal stability 0] . In the following, we consider 
the general case D > 4 from a unified point of view. 

A pole of the first type is related to the singularity of the function eq. I|34|) at w — 2 (the root, for details see the 
Appendix). We define the phase and amplitude via the integral l|31|) 

Y D (w,E) =I D (w,E)exp(i# D (w,E)). (40) 

It is not difficult to show that this pole emerges at arbitrarily small disorder for a negative phase i?d(2, E) < 0. 
The equation r i?r>(2, Eq) — defines the boundary of the region > Eq, where the physical solution is absent 
and, therefore, any disorder transforms the delocalized states into localized ones. For high-dimensional systems the 
delocalized states transform into states with non-exponential localization. The corresponding Eq values are Eq = 3.915 
(D=4), E = 4.365 (D=5) and E = 4.578 (D=6). 

For the region \E\ < Eq there exists the physical solution with a < a' (E), where <j'q{E) is a second threshold 
disorder value. The behaviour of this curve <j'q{E) (Fig. lb) is determined by resonance phenomena. The integral in 
(|31|l consists asymptotically of a power function, t^ ^ 1 ^ 2 , and a product of trigonometric functions. E.g. the function 
cos(2i — 7r/4) comes from every Bessel function. If we represent this product as a sum of monochromatic waves, then 
we denote the existence of a wave with zero frequency as a resonance. For w = the first non-trivial resonance lies 
either at E c = 2 (D=4,6,...) or at E c = 4 (D=5,7,...)- 

For \E\ < E c a pole of the second type appearing at higher levels of disorder a > cr' (E) must be considered. This 
type of pole emerges at purely imaginary values of the parameter w — iv. It corresponds to the roots of the equation 



a 2 Q-(v) = 1, 



(41) 
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FIG. 1: Threshold disorder values: (a) ao(E) for the transition from the non-exponential to the exponential localization and 
(b) a'o(E) for the transition from the delocalized to the (non-exponentially) localized states. The curves are enumerated with 
the values of D. 



where 

1 f 00 

Q-(v) = —= = / [J (2t)) D ~ 1 cos(Et)exp(-vt)dt. (42) 

V4 + v 2 Jo 

A physical solution is acceptable here only if the poles lie on the unit circle, \z\ = 1 or w E [0,2] for our sector w. 
A marginally stable solution corresponds to the value w — (or v — in eq. (|42|) '). The threshold disorder value is 
given by 

^(£)=Q_(0)- 1/2 . (43) 

For g < (t' q (E) the function H-(z) does not possess poles, here we find the region of stability of the extended states. 

Spaces of dimension D — 4, 5 possess a certain sensitivity with respect to the resonance phenomena mentioned 
above. For E c < \E\ < Eq the line of the poles exhibits a deviation from purely imaginary values of the parameter w 
and touches the real axis at the point u' G [0, 2]. This point corresponds to the condition for the phase of the integral 
*&d(u' ,E) — and also belongs to the unit circle (marginal stability). The corresponding threshold disorder value we 
denote again as a' (E): 

,„ {E) Ji^.\'\ (44) 



We see that the function a' (E) is in general singular in the energy. In going from the energy value E c to the value Eq 
the parameter v! increases monotonically and reaches finally the value u' — 2. The function cr (E) goes continuously 
to zero for \E\ — > E . 
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For dimensions D > 6 there are no resonance phenomena. These resonances are so weak that only one equation 
(|43H remains valid in the whole range of energies \E\ < Eq. I. e. although D c = 4 is a certain critical dimension (here 
non-exponential localization arises), a qualitative agreement of all results is only obrtained for D > 6. 

In this way it emerges from our exact analytic theory that for Anderson localization the upper critical dimension 
is Dq — 6, i.e. the problem in a certain way shows a similarity to percolation |30j - and this may not be totally 
unexpected. However, we also note that as a rough estimate for the upper critical dimension the value D c = 4 
may also be accepted. It corresponds to a first and important step of the qualitative saturation of the results (the 
possibility of existence of non-exponentially localized states), whereas a complete saturation only obtains at D=6. 

Quite generally in the range \E\ < E and under the condition a < a' (E) there exists a region of existence of 
stable delocalized states. Because in this range also non-localized states coexist with the other ones, this means that 
in this range the good metal (delocalized states) and the bad metal (non-exponentially localized states) can form a 
heterogeneous system, whose properties depend on the relative proportions of the subsystems; e.g. the bad metal, if 
the subsystem of non-exponentially localized states percolates. Fig. presents results of a numerical solution of the 
resulting equations. 

IV. LOW-DIMENSIONAL SYSTEMS 
A. Analytic solution 

Independently from whether D c — 4 or Dq = 6 are taken as the upper critical dimension, it is necessary to consider 
3-D systems as low dimensional ones. We approach here a widely held opinion that all states in D = 3 should be 
stable against perturbations and should have a finite radius of convergence. This, however, is not correct and would 
only be valid if the whole field of the scaling theory of localization would be faultless. 

The eci. (|37|l has always a solution for the physically important cases D = 2,3. We refer to these cases as 'low- 
dimensional' ones. For the low-dimensional systems the integral eq. (|38|) can be evaluated analytically (the formulas 
can be found in the tables of Laplace transforms). The corresponding pole diagrams are shown graphically in the 
Appendix. The system function has a pole at z = A — exp(27) with 7 > 70, 2 cosh(7o) = uq = 2p+ \E\ and p = D — 1. 
Note, however, that the above mentioned feature of the low-dimensional systems is caused by the divergence of the 
integral f2+(ito). 

If one applies the obtained equations 1)43(1 for <tq(E) to low-dimensional systems with D — 2,3, then formally 
<7q{E) = 0. Therefore, even infinitesimal disorder leads to solutions with exponential localization. Because the curve 
<Jq(E) forms the border between exponentially localized and non-exponentially localized states, this simply means that 
in this model for low-dimensional systems non-exponentially localized states are impossible (see, however, a discussion 
below). 

For such systems the second curve Oq{E) takes over the role of the mobility edge. The shape of this curve (Fig.2.b) 
for low-dimensional systems has a certain similarity with the ones for high-dimensional systems. Here we must stress 
that cf'q(E) does not represent a true mobility edge, but an upper limit for the coexistence of the two phases. 

For D = 2 in 8] it has been found analytically that a' (E) = 2(1 - £ 2 / 4 ) 1/4 for E < E o where E a = 2 = D. I.e. 
in this case the value Eq corresponds exactly to the midst of the band half-width 2D = 4. For E > Eq any disorder 
transforms the delocalized states into localized ones. 

For D = 3 the value Eq = 3 = D corresponds again to the midst of the band half- width 2D = 6. The resonance 
value mentioned above E c lies for D = 3 at E c = 4. Because Eq < E c , the resonance does not play a role here. The 
mobility edge <j' (E) again follows from eq. (|43|) . 

The case of the 3-D system shows an exception in the energy range Eq < \E\ < E' = 3.367, where the function 
(t'q(E) exhibits so-called reentrant behaviour. The proof of this requires a special investigation (see pole diagrams in 
the Appendix). 

Note that for the 3-D system in the band center cr (0) = 0, because the integral in eq. diverges. Fig. ^presents 
results of a numerical solution of the resulting equations. These results are compared to the analytical result of @ 
for the two-dimensional system (curve 2). 

B. Phase diagram and logarithmic divergence 

Although the results for D = 3 are formally exact, they require a special discussion. The mathematically correct 
results are physically acceptable only if they are stable against small perturbations (e.g. small changes in the model 
definition, fluctuations of parameters in the model). From this point of view two results have to be questioned: (i) the 
singularity of cr' (E) in the band center, <7q(0) = 0, and (ii) the non-existence of non-localized states, ctq(E) = 0. Point 
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FIG. 2: Threshold disorder values: (a) ao(E) for the transition from the non-exponential to the exponential localization and 
(b) cr' (E) for the transition from the delocalized to the localized states. The curves are enumerated with the values of D. The 
curves 3' and 3" correspond to the model eg. l45H for k = 0.01 and K = 0.1. 



(i) arises mathematically as a consequence of the logarithmic divergence of the integral 1(42(1 for v = and \E\ — ► 0. 
Point (ii) arises again via a logarithmic divergence of the integral 1(38(1 for u = uq. I.e. basically, the case D = 3 is 
nothing else but a type of logarithmic deviation from the high-dimensional case. Under certain conditions one can 
find the perturbations which are capable of regularizing the mentioned logarithmic diver gen ce, i.e. to transform them 
into a finite term. One could e.g. surmise that in this regularization correlated disorder [3l| might play an important 
role. One must stress here that we are dealing only with results for tight-binding Hamiltonians with diagonal disorder. 
It is largely unclear, whether this property remains valid also for non-diagonal disorder. 

To illustrate this let us first consider a purely mathematical problem. Is it at all possible to confirm the results (i) 
and (ii) cither by different numerical computations or analytical approaches? The answer is - no! Any deviation from 
the exact solution (via numerical or analytical ways) automatically generates results of a mean-field theory. A mean- 
field theory has a certain qualitative agreement with the results of the exact theory, but only for high-dimensional 
systems. Everything looks as if any uncontrolled deviation from the exact theory has added additional dimensions to 
a 3-D system. As a mathematical model let us further consider a 4-D system, where, however, a coupling involving 
this additional dimension is exceedingly weak (parameter /t<l). We start from a generalization of the function 120(1 
for D = 3: 

2 

£(k) =.E-2^cos(fcj)-2Kcos(fc 3 ). (45) 

3=1 

Here the term involving k corresponds to the hopping matrix element into the fourth dimension. Figj^J gives the 
numerical results for 2 cases: k — 0.01 (curve 3') and k = 0.1 (curve 3"). One clearly sees that all results, which have 
nothing to do with the logarithmic divergence, e.g. the entire curve cr' (E) with the exception of the point E = 0, 
remain extremely stable. Even the change in the parameter Eq is of the order of 0(n 2 ). 
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The 'logarithmic' results on the other hand turn out to be completely unstable. Even the small values of the coupling 
(or regularization) parameter n produce results which are in qualitative agreement with those for high-dimensional 
systems. The curves 3' and 3" look as if they were an extrapolation of the corresponding curves for (Jq(E) from 
Fig- ffli. Here the inequality <tq(E) > cr' (E) also applies. Reducing the parameter k slowly moves the mobility edge 
<Jo(E) downwards, because the dependence on the parameter n is extremely week (logarithmic), cr (E) ~ l/ln(K _1 ). 
Only for even smaller values of the parameter k <C 0.01 one might perhaps see traces of the exact results, because in 
this case one has ao{E) < a' (E) for \E\ < E - 

V. CONCLUSION 

Although the formal investigation of the Anderson model of localization (tight-binding Hamiltonian with diagonal 
disorder) for higher spatial dimensions D might at first look very abstract, the exact analytical results supply us with 
clear physical consequences. The Anderson problem of localization and the percolation problem belong to the same 
class of critical phenomena: both have the same lower and upper critical dimensions. I.e., although the Anderson 
model appears to be much more complex and richer, certain fundamental results appear to be transferable. Percolation 
is possible for 2-D systems, this corresponds to the existence of a metal-insulator transition in disordered 2-D systems. 
In this sense there is no reason to believe that the existing contradictions between theory and experiment for 2-D 
systems point to an incompleteness of the Anderson model. On the contrary, our analytical investigation has shown 
that a tight-binding Hamiltonian is presumably sufficient for this purpose. 

The main problem of the theory thus does not rest in the Hamiltonian, but rather in the interpretation of the results, 
which mainly derive from numerical work. Our analytical and exact results demonstrate the necessity of interpreting 
the phase transition in the framework of first order phase transition theory and this holds independently of the spatial 
dimension D > 2. If one, however, attempts to apply procedures which have only been developed for systems with 
a second order phase transition (and this is the general case), one does not necessarily obtain wrong numbers, but 
an incomplete or even wrong interpretation. See an example in 8], where an analytical (exact) scaling function for 
the 2-D system has the same form as obtained by numerical scaling. The physical interpretation is, however totally 
different. We hope that the necessary corrections of the numerical tools are possible to detect the first order phase 
transition. 

The rather large value of the upper critical dimension for the Anderson localization (and the percolation) problem 
permits to consider 2-D and 3-D cases as low-dimensional systems. Thus a revisiting of the results is also necessary 
for 3-D systems. We have found that the 3-D case is nothing else but a type of logarithmic deviation from the 
high-dimensional case. As a consequence results a certain instability of the results, whose details are discussed in the 
text. 

We give here a short summary of the main results. 

• For the Anderson localization problem there exists an upper critical dimension Do = 6. This value is also 
characteristic of a related problem: percolation J2ij. For D > Dq all phase diagrams are qualitatively the same, 
only the corresponding critical values develop in a monotonic way. One can also say that this is the property of 
a mean-field theory, although in this case a mean-field theory does not exist as a closed theory. 

• There exists also a second upper critical dimension D c — 4, which has a different meaning. The states with 
non-exponential localization are formed only for D > D c , whereas for D < D c localized states are always 
exponentially localized. This second upper critical dimension D c divides the dimensions into two classes: high 
dimensions with D > D c and low dimensions with D < D c . 

The results for the nontrivial spatial dimensions D > 1 can be summarized as follows. 

(i) All states with energies \E\ > Eq are localized at arbitrarily weak disorder. The value of Eo depends on D and 
lies inside the band Eq < 2D. 

(ii) For \E\ < Eo states are only localized if the disorder a exceeds a critical value a' (E), otherwise a two-phase 
system is formed from an insulating and a metallic one. This differs from the traditional point of view which considers 
the localization transition as a continuous (second order) transition. Should the standard interpretation of this system 
in the framework of first-order phase transition theory be applicable (which still has to be investigated) one can expect 
that qualitatively it has similar properties as other heterogeneous two-phase systems (e.g. the coexistence of water 
and ice). Then percolation problems might be important. 

(iii) Exponential localization always exists for the physically important cases D — 2,3. Non-exponential localization 
occurs only for higher dimensions D > D c = 4. In this case for \E\ < Eq and a < cr Q (E) first the heterogeneous 
system appears, where the difference between the two phases may be small (this has to be investigated). For <tq(E) > 
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a > ct'q(E) one finds a homogeneous system with non-exponential localization. Only for a > <7q(E) a system with 
exponential localization appears. 

(iv) ct'q(E) is in general not an analytic function of the energy E. There exist certain resonances. 
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APPENDIX A: POLE DIAGRAMS 

Parametric representation of the pole diagram 

A filter h n is characterized by a pole diagram of its image H(z). The principal definition ea. l|24|) . together with the 
two other ones, eqs. H30|) .H33 [1 . supply us with the parametric w-representation of the pole diagram, ea. 1)34(1 . Let us 
rewrite this relation into the form 

-L^ = l-a 2 R(w,E), (Al) 

where R(w, E) is generally a complex function of the complex variable w and energy E. The main idea is quite simple. 
The function H{z) has its poles where 

a 2 R(w,E) = 1. (A2) 

An elementary requirement for this is the condition on the argument of the complex function 

axgR(w,E) = 0, (A3) 

because a is positive. For the given energy value E ea. (|A3|) defines one or more lines (pole lines) in the complex w- 
plain. As we have already discussed in the text, for reasons of symmetry it suffices to analyse only a sector w — u + iv 
with u > and v > 0. Let us analyse one of these pole lines. For each point w on this line the position of the pole is 
determined via ea. (|32[l . The corresponding value of the disorder a is found from ea. l|A2() : 

a = R{w,E)- 1 ' 2 . (A4) 

Because for each value of w on the pole line there is associated a value of a, it is possible to indicate by an arrow next 
to the line in the diagram in which direction the poles move with increasing disorder. 

Filter H+(z) 

Let us consider first the simplest case, the filter H{z) = H + (z). Here the poles are connected with the notion of 
the Lyapunov exponent. In the parametric representation there exists a simple relation 

w = 2cosh(7); (A5) 

thus the eas. i|A4|) . l|A5|) together supply a connection between a and 7. 

For a 2-D system it is possible to analytically evaluate the corresponding function R(w, E) @: 

R(w,E) = - 7 ^=[ 1 ^ 2 — + n = =fL (A6) 

2\/w 2 - 4 y/(w + E) 2 - 4 y/(w - E) 2 - 4 



or 



R(w,E) = - ? 2==x (A7) 
2Vw 2 - 4 

1 



'y/(w + E + 2)(w + E-2) yj(w- E + 2)(w- E-2) 
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We clearly see that the function (|A7|) possesses for the given energy several points, where it diverges. These values 
w = — 2 — E, w = 2 — E, w = E — 2 and w — E + 2 are real. They are nothing else but the resonances discussed 
in the text. To detect these we investigate the w- parametric representation of the integral ea. (|31[l . Every Bessel 
function Jo(2i) contributes asymptotically a trigonometric function, cos(2i — tt/4). In addition there exists another 
energy dependent trigonometric function, cos(Et). If we represent all these functions via complex exponentials, we 
obtain under the integral in eq. (|31|1 asymptotically a product of the power function and a sum (with well determined 
coefficients) of exponents exp[i(w — Wj)t\. These are the resonance values uij. In the 2-D case we have a strong 
resonance (the function R(wj,E) diverges). 

It is easy to establish that the function R(w,E) satisfies eq. (|A3(l for the real value of the parameter w under the 
condition 

w > uq = max{wj}, (A8) 

where for the 2-D system uq = 2 + \E\. I.e. here exists a line of poles which starts from the point uo = 2 + \E\ (it 
corresponds to the value a = 0, because here R(uq, E)^ 1 / 2 = 0) and moves on with increasing value of the parameter 
a always along the real axis. It is also easy to find that in this case no other pole lines exist. This result can be 
generalized for higher dimensions. So one can establish that the value of uq found from ea . 1)36(1 corresponds precisely 
to the condition ea. (|A8fl : we always have to deal with the same resonance. 
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FIG. 3: Parametric representation of the pole diagram for an unstable filter: a) low-dimensional systems with D = 2, 3; b) 
high-dimensional systems with D > 4. 



The structure of the pole diagram is nearly identical for all dimensions D, Fig03 There exists only one pole line, 
which originates at the point w = uq. The direction of increasing a is denoted by an arrow. The only difference is that 
for low-dimensional systems the starting point uq corresponds exactly to a = (this is marked in the diagram with 
a black circle), for high-dimensional systems the value (Jq{E) = R(uq, E)^ 1 / 2 is finite (the point w = uq is marked 
with a white circle). In the latter case it is impossible to find for smaller values of the disorder a < ao(E) a point in 
the diagram, where the condition ea. ((A2(l is fulfilled. The filter H+(z) possesses no poles. The border thus found via 
ea. ()39|l has been defined as mobility edge. 



Filter J2_(z) 

This filter has a physical interpretation only under the condition Q that the function H-(z) has either no poles or 
that they belong to the unit circle, \z\ = 1. In the parametric w- representation the unit circle (in the sector w = u + iv 
with u > 0, v > 0) corresponds to an interval on the real axis u £ [0,2]. I.e. it is necessary first to find the single 
points or even lines in this interval u € [0,2], where the condition ea. (|A3|l holds. Later on one must also find out, in 
which way the pole line leaves this interval. 

It is easy to ascertain that one of the possible pole lines always lies on the imaginary axis v, and in the interval 
v > 0. On this line (trivial pole line) only the point v — could have a physical interpretation, because the point, 
w = 0, also belongs to the unit circle. The existence of other pole lines strongly depends on the space dimension D 
and the energy E. 

As the simplest example let us first consider the filter for the 2-D system , where the corresponding function has 
an analytical form: 

R(w, E) = 1 [ 1 + / ]. (A9) 

2V4 - w 2 ^4 - (w + E) 2 y/4-(w- E) 2 
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Because the interval u € [0, 2] plays an important role for the physical interpretation, we consider first the resonance 
values {wj}, which lie precisely in this region. It is easy to calculate that for each energy value only one resonance 
w = uq exists. For energies < \E\ < 2 we have Uo = 2 — \E\. In the remaining energy region, 2 < \E\ < 4, one has 
uq = \E\ — 2. In the 2-D system the resonances are so strong (the function R(w,E) diverges at w = uq) that these 
determine in a unique manner the beginning of the new pole lines and their direction, Fig0] 




FIG. 4: Parametric representation of the pole diagram for a 2-D system: a) energy range < \E\ < 2; b) energy range 
2 < \E\ < 4. 



In the interval < \E\ < 2 the pole line starts exactly at the point uo = 2 — \E\. This point corresponds to the 
value of the disorder a = and is denoted by a black circle in the diagram. With increasing value of a the pole line 
follows exactly the real axis u, until it reaches the point u = (w = 0). At this point the corresponding disorder 
i?(0,i?)~ 1 / 2 is finite (white circle in the figure), and we have denoted this in ea. (|43|l as a' (E). Because this line 
is always a unit circle (marginal stability), we have physical solutions. These exist only under the condition that 
< g < a' (E). If the disorder crosses the border a' (E), the pole line leaves the point w — and follows further 
the trivial pole line (imaginary axis). Here, however, a physical interpretation is no longer possible and consequently 
there are no extended states in the range a > <r' (E). 

In the interval 2 < \E\ < 4 the behaviour of the pole line determines the resonance at w = uo = \E\ — 2. This line 
again starts at the value uq, which corresponds to the disorder a = (black circle in the figure). Then this line leaves 
the real axis and takes its course into the complex plain until it reaches a particular point vq (bifurcation point) on 
the imaginary axis v. The later path makes use of branches of trivial pole lines. The directions belonging thereto 
are again marked with an arrow. Thus we have for a > no points in the phase diagram which can be interpreted 
physically. I.e. an infinitesimal disorder already suffices in this energy range to destroy all extended states. One may 
also write that in this range (t' (E) = 0. 

Formally between the 3-D case (low-dimensional system) and the cases D > 4 (high-dimensional systems) there is 
a quantitative difference. For high-dimensional systems one does not find a divergence of the function R(w, E) at the 
resonance w — Uo, but only a jump in the argument of this complex function. For the 3-D system this divergence 
exists, on the other hand, but in contrast to the 2-D system the resonance is very weak (logarithmic divergence). 
In all these cases the resonances have no direct influence: they do not generate a pole line emerging from the point 
uq. There is, however, an indirect influence of the resonances, because every resonance defines the argument of the 
function R(w, E). 

In fig'03 different types of pole diagrams are presented. Fig|5^ is typical for small values of the energy. Although 
the resonance at w — uq exists (this point is marked in the figure with a black square), the weak resonance does 
not generate a new pole line. There remains only the trivial pole line, where only the point w — has physical 
significance. The value R(0, E)~ x l 2 defines again the function cr (E) already mentioned. Physical solutions exist 
only for a < <Jq{E) (no poles on the unit circle). FigOJ) is typical for large values of the energy. In the parametric 
representation of the filter function, ea. (|34|l . one can in addition see a point, where the function R(w,E) diverges. 
This is the point w = 2 (a type of energy independent resonance, which corresponds to the factor V4 — w 2 ). One can 
see that this point takes over the role of the resonance at uq. The pole line starts from the point w = 2 (this point 
corresponds to the value a = and is therefore again marked with a black circle). The continuation is similar to 
FigEJ): the line continues to the imaginary axis and reaches it at the point v = vq. Although the points w — uq and 
w = 2 are not identical, FigQJ} and FigOj are qualitatively similar, and they have the same physical interpretation. 
In this range of the energy an infinitesimal disorder destroys all extended states. 

An indirect influence by the resonances consists in the fact that one can change in their neighbourhood the argument 
of the complex function R(w, E) in such a way that the condition ea. !A3(l becomes valid at a point w — u\ see FigEt 
and FigOJl. The point u' is not a resonance, here the function R(u', E) remains finite (white circle in the pole 
diagram) . 
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FIG. 5: Different types of pole diagrams for the systems with spatial dimensions D > 3. For details see the text. 



Let us consider further on two examples. In the 3-D case one can divide the energy values into three ranges. In 
the range < \E\ < 2 the resonance is found at uq = \E\, which, however, has no importance. We have the pole 
diagram of Fig[S^. In the range 2 < \E\ < A the resonance occurs at uq = 4 — \E\. Here the above mentioned point v! 
arises, but under the condition \E\ > Eq = 3.00. I.e. in the range 2 < \E\ < Eq the pole diagram of Fig|3^, remains 
valid. In the range Eq < \E\ < 4 we find a different typ of diagram, Fig[!J;. The point v! corresponds to the condition 
< v! < uq. For E — ► 4 the point u' moves towards the value v! = 0, and thus arises the bridge between the new 
pole line and the trivial pole line. After this the diagram is qualitatively the same as FiglSJ}: there are no extended 
states. This is also valid in the range 4 < \E\ < 6, where uq = \E\ — 4. 

Figure describes a complicated case which does not have an unambiguous interpretation. Formally this is the 
only diagram which has two pole lines. The other ones consisted always of a single pole line, although bifurcation 
points were also possible. For infinitesimally small disorder the resonance at w = 2 is important, as in the diagram 
of FigUJ). Here with increasing disorder a, the pole line which starts at the point w = 2 (black circle), leaves the 
interval u £ [0, 2). The corresponding poles have no physical interpretation, which corresponds to the annihilation of 
extended states via infinitesimal disorder. If the disorder increases further this pole line approaches again the interval 
u e [0, 2], and reaches it at the point v! (Fig. |SJ:). The corresponding value of the disorder o\ = R(u' , E)~ x l 2 is finite 
(white circle). The physical interpretation now depends to which value of the disorder 02 — -R(0, E)~ x / 2 corresponds 
the point w = (white circle) on the trivial pole line. 

If o 1 < 02, then a sort of gap arises in the disorder in such a manner that in the range <r\ < a < 02 no values of 
the parameter w correspond to the pole line. I.e. the filter H(z) has no poles in this range. Consequently a physical 
interpretation of the solution is possible here and we obtain the reappearance of extended states at finite disorder 
values in this special interval. This condition, o\ < 02, is valid for 3-D systems only in a very narrow range of the 
energy, Eo < \E\ < E' , where E' — 3.367. The behaviour of the curve (j' (E) is shown in Fig|5] (so-called reentrant 
behaviour). We can clearly see that to each energy value in the range Eq < \E\ < E' Q are associated three values 
of (t'q(E); these are in particular cr' (E) = - localization via infinitesimal disorder, cr' (E) = o\ - reappearance of 
extended states, and cr' (E) — 02 - again localization. 

Outside this range, E' < \E\ < 4, one has u\ > 02. For this condition there are no points on the pole line which 
permit a physical interpretation. I.e. after this annihilation of extended states via infinitesimal disorder it is not 
possible for extended states to reappear. 

In the 4-D case we find a different sequence of resonances. In the range < \E\ < 2 the resonance occurs at 
uq = 2 — \E\, and the diagram Fig[5k is valid. In the range 2 < \E\ < 4 the resonance is found at uq — \E\ — 2. Here 
the point u' (uo > u' > 2) arises. The pole diagram corresponds to[3l. The pole line starts from a white circle, which 
corresponds to the value cr' (E) > 0. Only this point in the pole diagram has a physical interpretation, and it means 
the existence of extended states under the condition of small values of the disorder parameter, a < <r' (E). For E — > 4 
the values u' and uq move towards w = 2. Because always u' > uo, this means that the point v! reaches the value 
w = 2 before the point uq, and this occurs for Eq = 3.915 < 4. For E > Eq the diagram |Sj3 is still valid, although 
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the series of resonances is not yet exhausted: a resonance arises at uq = 6 — \E\ for 4 < \E\ < 6, and another one at 
i*o = \E\ — 6 for 6 < \E\ < 8. Here also extended states are not possible. 
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